Anderson localisation in tight-binding models with flat bands 
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We consider the effect of weak disorder on eigenstates in a special class of tight-binding models. 
Models in this class have short-range hopping on periodic lattices; their defining feature is that 
the clean systems have some energy bands that are dispersionless throughout the Brillouin zone. 
We show that states derived from these flat bands are generically critical in the presence of weak 
disorder, being neither Anderson localised nor spatially extended. Further, we establish a mapping 
between this localisation problem and the one of resonances in random impedance networks, which 
previous work has suggested are also critical. Our conclusions are illustrated using numerical results 
for a two-dimensional lattice, known as the square lattice with crossings or the planar pyrochlore 
lattice. 
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I. INTRODUCTION 

Hopping problems on certain frustrated lattices have 
the unusual feature that some of their energy bands are 
dispersionless. Systems of this kind have been studied in 
several different contexts, including both itinerant ferro- 
magnetism within the Hubbard models and frustrated 
antiferromagnetism within the Heisenberg modeli^ In 
the first of these examples, the divergent density of states 
for a flat band stabilises ferromagnetic order by ensuring 
that the Stoner criterion is satisfied. In the second exam- 
ple, the existence of a flat band in the hopping problem 
is closely linked to macroscopic ground state degeneracy 
in a classical antiferromagnct defined on the same lattice. 

The degeneracy of states in dispersionless bands is eas- 
ily lifted by perturbations. In this paper we examine the 
effect of one such perturbation: a weak random poten- 
tial. For a flat band of the clean lattice, weak disorder 
sets the only energy scale and controls the localisation 
properties of eigenstates. Previous work in this area^ 
while identifying the weak-disorder limit as interesting, 
has focussed mainly on Anderson transitions occuring at 
finite disorder strength. We expect transitions occurring 
at finite disorder strength to be in the same universal- 
ity class as those on unfrustrated lattices, and anticipate 
that special features of fiat-band localisation occur only 
in the weak-disorder limit. 

The weak-disorder flat-band localisation problem 
shares an obvious feature with the one of Anderson lo- 
calisation in integer quantum Hall systems, in the sense 
that a Landau level can be viewed as a continuum ana- 
logue of a flat band. For a Landau level, a property 
of the projection operator onto the level - its non-zero 
Chcrn number - ensures that disorder cannot localise all 
states^ Remarkably, we find for frustrated hopping prob- 
lems that localisation is also controlled by a feature of 
the projection operator onto the flat band. In the cases 
we are concerned with, its matrix elements in real space 



decay as a power law of separation, and this power is 
equal to the spatial dimension. As a consequence, a ran- 
dom potential projected onto the flat band gives rise to 
a localisation problem similar to one in which the mag- 
nitude of hopping matrix elements between distant sites 
decreases as a power of their separation. A system of this 
kind is known to be critical when the power is the same 
as the spatial dimensional^ and the one-dimensional case 
of power-law banded random matrices has been studied 
extensively^ On these grounds, we expect eigenstates in 
frustrated tight binding models with weak disorder to be 
critical, and we use a numerical study of their multifrac- 
tal properties to demonstrate this for a two-dimensional 
example. Separately, we show that the flat-band locali- 
sation problem can be mapped onto the problem of res- 
onances occurring in random impedance networks, for 
which cigenfunctions have also been found to be critical 
in an earlier numerical study^ 

While our exclusive focus here is on flat-band localisa- 
tion, this problem has close links to a number of other 
topics of current interest. The most direct connection 
is to a geometrically frustrated antiferromagnet that has 
its spins polarised by a strong applied field: its single 
spin-flip excitations are described by a frustrated hop- 
ping Hamiltonian ; 12 i 13 and in the presence of weak dis- 
order may exhibit the features we discuss in this paper 
(although the form of disorder we treat is site-diagonal, 
while randomness in the strength of exchange interac- 
tions would appear as off-diagonal disorder). More gen- 
erally, the mathematics of flat bands is closely linked to 
the physics of Coulomb phases^ in classical statistical 
mechanical systems - notably dimer models^ and frus- 
trated magnets^ There has recently been a general ef- 
fort to understand transitions out of the Coulomb phase, 
induced by perturbations^^— In particular, flat-band 
localisation can be viewed as a linearised version of the 
transition induced by weak disorder in a frustrated mag- 
net, between the Coulomb phase and spin glass order 2& 



Finally, we note recent studies of Bose condensation in 
flat bands^ and of the combined effects of disorder and 
interactions in the Falicov-Kimball model on a lattice 
with a flat band^ 

In Section [TT] we define the models we consider, dis- 
cuss projection onto the flat band and set out a mapping 
between the fiat-band localisation problem and random 
impedance networks, while in Section IIIII we present re- 
sults from simulations. 



II. MODELS AND THEORETICAL 
BACKGROUND 

Several parallel terminologies have been used to de- 
scribe the tight-binding models we study here, which 
all involve a single non-zero hopping energy t between 
nearby sites labelled i at positions rv on a periodic lat- 
tice in d dimensions. In informal terms, for each model 
the sites can be grouped into clusters, labelled a, /3 . . ., in 
such a way that: (i) each site is shared between two clus- 
ters; (ii) every pair of sites belonging to a given cluster is 
linked by hopping; and (Hi) there is no hopping between 
sites from different clusters. A second lattice appears 
naturally from this construction: it has its sites located 
at the centers of the clusters of the first lattice, and the 
sites of the first lattice lie at the mid-points of bonds of 
the second lattice. As described by Mielkei using the 
language of graph theory, if the second lattice is denoted 
by the graph G(V, E) with vertices V and edges E, then 
the first lattice is the line graph L(G) of G, which has 
as its vertex set Vl the edge set E of G, and has edges 
El between two vertices in Vl if and only if the corre- 
sponding edges in G share a vertex. Our study is hence 
of tight binding models on line graphs. Alternatively, G 
is known as the simplex lattice^ or parent lattice^ and 
L(G) is called its medial lattice. One distinction within 
this class is important to us: for reasons indicated below, 
we restrict ourselves to simplex lattices that are bipartite 
and expect different behaviour in the opposite case. We 
note that some rather more general classes of tight bind- 
ing models with flat bands have also been considered^ 
which we do not examine here. 

A two-dimensional example of the class of lattices 
which we are interested in is the square lattice with cross- 
ings, also called the planar pyrochlore lattice, shown in 
Fig. Q] For this the simplex lattice is a square lattice. A 
second two-dimensional example is the kagome lattice^ 
which has the hexagonal lattice as its simplex lattice. 
A three-dimensional case is the pyrochlore lattice^ with 
the diamond lattice as its simplex lattice; another is the 
octahedral lattice^ with the simple cubic lattice as its 
simplex lattice. 

When specifying the tight-binding Hamiltonian T-i^ 
for the clean system it is convenient to include a constant 
term, chosen so that the flat bands are located at energy 
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FIG. 1: The square lattice with crossings, or planar py- 
rochlore lattice. Filled circles denote lattice sites and lines 
represent hopping matrix elements of magnitude t in 1v \ 

zero. The Hamiltonian then has elements 

%(°) _/*[! + J e same cluster ,^ 

W 1 otherwise . ^ ' 

Its expectation value in a state with amplitudes ipi can 
conveniently be written in the form 

Mtt (0 V> = *EIZ>l 2 - (2) 

Taking t > 0, it is clear from Eq. ^ that the eigenvalues 
of are non-negative. It is also evident that states 
having 

for all a and {<p\tp) ^ are eigenstates with energy zero. 
Such states can be constructed on the lattices we study 
by taking ifi = ±1 with alternating signs on a closed 
loop of neighbouring sites that includes exactly two sites 
from every cluster visited by the loop, and setting tp^ = 
at sites not lying on the loop. (The construction with 
alternating signs requires the closed loop to have even 
length. Since the loop can be viewed as a walk on the 
simplex lattice, this property is guaranteed provided the 
simplex lattice is bipartite, but not otherwise.) These 
loops can be chosen to be finite and short: for example, as 
squares on the planar pyrochlore lattice or as hexagons on 
the kagome lattice. Translated copies of such eigenstates 
form a band of zero-energy states, which is the flat band 
we are concerned with. 

In what follows, properties of the projection operator 
V onto the zero-energy states play an important role. 
When a band is separated by energy gaps from other 
bands, the real-space matrix elements Vij decay expo- 
nentially at long distances as a function of separation 
|rj — rj| between sites^ However, for models in the class 
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we have defined (here it is again crucial that the simplex 
lattice be bipartite), the zero-energy band touches a dis- 
persive band at one or more points in the Brillouin zone^l 
and Vij decays more slowly with separation. For any spe- 
cific lattice its form can be computed explicitly from the 
eigenvectors of H^°\ but to see the generic long-distance 
properties it is better to move to a continuum treatment. 
To this end, following discussions of the Coulomb phase 
in dimer models^ and antiferromagnets^ we first intro- 
duce real space unit vectors aligned along each bond i 
of the simplex lattice, and all directed from one sublat- 
tice (chosen as a matter of convention) towards the other. 
Then we associate a (i-component vector field B(r^) with 
a set of amplitudes tpi, by defining 

B(r;) = tpiBi . (4) 

The zero energy condition, Eq. (|3|), is also the condi- 
tion that the field B(r.;) has zero lattice divergence at 
every site of the simplex lattice. To understand the long- 
distance behaviour of V, we therefore consider the projec- 
tion operator onto divergencelcss fields that are functions 
of a continuous position variable r. The reciprocal space 
form of V is block-diagonal in wavevector q. On the lat- 
tice the blocks act on the space of sites within a unit cell, 
while in the continuum they are d x d matrices acting on 
the components of B, with the form 

q 2 Si m - qiq m 

— ? — ■ (5) 

In real space this yields 

„ dnr m ~r 2 5i m 

l J R,R+r - C d , (6) 

where Cd is a dimension-dependent constant. As ad- 
vertised, the projection operator falls off as r~ d in d- 
dimensions. 

Returning to the tight-binding model, we introduce 
disorder as a site-diagonal random potential with matrix 
elements 

«g } = Sis* ■ (?) 

The Vi are independent, identically distributed random 
variables having zero mean and unit variance. Combining 
ingredients, the localisation problem we study has the 
Hamiltonian 

u = n w + n w (8) 

and we arc interested in the weak disorder limit, obtained 
by taking t —> oo. 

Consider for T~L the density of states in energy E when 
t > 1. At \E\ > 1 it is similar to that of H {0) , but for 
\E\ < 0(1) it is dominated by the contribution from the 
disorder-broadened flat band, which has energy width 
0(1). We wish to understand simplifications at large 
t in the eigenvectors of T~L belonging to this band. Let 



If) = \f\\) + IVJ-) be such an eigenvector, which we have 
separated into components that obey V\ip\\) = \<p\\) and 
(1 — V)\(p±) = \<p±) and lie respectively within and per- 
pendicular to the space spanned by the flat bands of the 
clean system. Because \tp\\ ) = 0, the eigenvalue equa- 
tion E\ip) = W\tp) implies 

(E-UW)\<p}=UW\tp ± ). (9) 

Since we expect \(E — 'H < - 1 " , )\'p)\ ~ 0(1), we require 
\\ip±)\ ~ Oir 1 ) to ensure that \H {0) \<p±)\ ~ 0(1) when 
t is large. At leading order, Eq. © therefore simplifies 
to 

{E-H^)\ n )=U^\<p±). (10) 

To make a connection with random impedance net- 
works, note that V and (1 — V) effect a lattice version 
of the Helmholtz decomposition: V projects onto lattice 
fields that are solenoidal, while 1 — V projects onto lattice 
fields that are irrotational. In consequence, 'H^\tp±_} can 
be expressed as the lattice gradient of a potential defined 
at the sites of the simplex lattice. Let a and (3 be the 
simplex lattice sites linked by the simplex lattice bond i, 
with e, directed from a to j3. Then there exist simplex 
site potentials $q, such that for every i 

<?:|H (0) |<^} = * Q (li) 

In addition 

<t\(E-HW)\ n ) = {E- Vi ){i\ n ) (12) 
and so with w a p{E) = (E — v.;) -1 we have 

(i\ n )=w a p(E)($ a -$p). (13) 
Now Eq. implies for every a that 

Y^v> a p(E)(*a = (14) 

p 

where the sum on /3 is over simplex lattice sites that are 
nearest neighbours to a. 

Eq. (fT4")l is also the equation that describes resonances 
in an impedance network^ It has non-trivial solutions 
for <E> Q only if E is an eigenvalue of H. In the impedance 
network interpretation, it is Kirchoff's law for current 
conservation at node a of the network, 4> a is the volt- 
age at that node, and w a p{E) is the complex conduc- 
tance, or inverse impedance, between nodes a and (3. 
This impedance is frequency-dependent and the equa- 
tion has non-trivial solutions at the resonant frequen- 
cies, which are real if the network is loss- free. Such 
resonances have been studied for a binary distribution 
of impedances on the square lattice by Jonckhcerc and 
LuckAi In particular, these authors find from numerical 
calculations that the electric fields associated with res- 
onant states, which correspond in our notation to the 
combination [E — Vi){i\ip\\), have a multifractal distri- 
bution. Their conclusion is entirely consistent with our 
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FIG. 2: The probability distribution of nearest-neighbour 
level spacings P(s), obtained for system sizes L = 34 (x), 
44 (*) and 54 (□). Inset: detailed view of behaviour at small 
s. The Poisson distribution is represented by the full line, and 
the Wigner-Dyson distribution by the dashed line. 



view of fiat-band localisation, as being automatically crit- 
ical because V is long range. Moreover, this connec- 
tion to flat-band localisation provides an explanation of 
why resonant states should generically be critical for an 
impedance network in any number of dimensions. 



III. NUMERICAL RESULTS 

To illustrate these ideas and study the approach to 
the large t limit, we have calculated eigenfunctions and 
eigenvectors for % on the planar pyrochlore lattice of 
Fig. [TJ We take the site potentials Vi from a Gaussian 
distribution and examine values of t up to 3 x 10 9 . We 
diagonalise H for square systems with periodic boundary 
conditions and side length L in the range 34 < L < 80, 
in units of the nearest neighbour spacing. 

The probability distribution P(s) of spacings s be- 
tween adjacent eigenvalues, measured in units of the 
mean spacing, is a standard diagnostic for localisation 
transitionSi 28 i 29 Metallic and localised phases are iden- 
tified by Wigner-Dyson and Poisson distributions^ re- 
spectively, while the critical point is characterised by a 
distinct, universal distribution^ We present in Fig.[2]our 
results for P(s). These were obtained using states from 
the central half of the flat band, and applying standard 
band-unfolding methods to compensate for the energy 
dependence of the average density of states. The data 
shown were calculated using t — 3 x 10 9 but results are 
essentially the same for any t > I0 3 . Within the sta- 
tistical precision of the data, the behaviour of P(s) is 
independent of system size. It is apparent from this fig- 



ure that the form of P(s) is intermediate between those 
for Wigner-Dyson and Poisson distributions, as expected 
in a system at a critical point. 

To probe the nature of eigenfunctions we have exam- 
ined the scaling with system size of averaged moments of 
the probability density. More specifically, cigenfunction 
fluctuations can be characterised by a set of generalised 
fractal dimensions 3 -^ if disorder-averaged moments have 
the power-law dependence 



|2ql 



= A q L~ T ^ 



(15) 



on system size, where r(q) is related to the fractal di- 
mension by r(q) = (q — l)D q and A q is a constant. For 
uniformly extended states D q = d and for localised states 
D q = 0, but for critical states a non-trivial dependence of 
D q on q is expected. As is standard^ we analyse multi- 
fractality via the spectrum of singularity strengths f(a), 
which is related to r(q) by a Legendre transformation: 
f(a) = qa — r(q) with q = d a f(a). It has the inter- 
pretation that L-^ Q ) is the measure of the set of points 
i at which the probability density in an eigenstate scales 
with system size as \ fi\ 2 ~ L~ a . We determine f(a) us- 
ing the procedure described in Ref. l32l . The results are 
shown in Fig. [31 The upper panel of this figure displays 
the size-depcndcncc of our estimates for /(a) at finite 
L, showing convergence to the large-L limit. The lower 
panel demonstrates that the form of f(a) that results at 
large L is t-independent provided t is large^ These final, 
large-L, large-i data are well- fitted by the parabola 



/(«) 



(a - ap) 2 
4(a - d) 



(16) 



with d = 2 and ao — 2.37. Correspondingly, we find 
£>2 = 1.27. Fractal properties in generalized random ma- 
trix ensembles depend on details of the model 3 ^ (for ex- 
ample, on disorder strength in power law random banded 
matrices^). Similarly, we expect fractal properties in 
weak-disorder flat-band localisation to vary with the 
choice of distribution for Vi, but we have not explored 
that aspect. 

In summary, these numerical results support our ex- 
pectation that weak disorder in a tight binding model 
with flat bands should give rise to cigenstates that are 
critical, in the sense of the Anderson localisation transi- 
tion. 
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